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Propagation of Ultrahigh Energy Nuclei in the Magnetic Field of our Galaxy 
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Abstract: In this work, we present detailed simulations for propagation of ultra-high energy (UHE) heavy nuclei, with 
E > 60EeV, within recent Galactic Magnetic Field (GMF) models [1,2]. We investigate the impacts of the regular and 
turbulent components of the GMF. We show that with UHE heavy nuclei, there is no one-to-one correspondence between 
the arrival directions of cosmic rays (CR) measured at Earth and the direction of their extragalactic sources. Sources can 
have several distorted images on the sky. We compute images of galaxy clusters and of the supergalactic plane in recent 
GMF models and show the challenges, and possibilities, of "UHECR astronomy" with heavy nuclei. Finally, we present 
a quantitative study of the impact of the GMF on the (de-)magnification of source fluxes, due to magnetic lensing effects. 
We find that for 60 EeV iron nuclei, sources located in up to about one fifth of the sky would have their fluxes so strongly 
demagnified that they would not be detectable at Earth, even by the next generation of UHECR experiments. 
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1 Introduction 

The Ultra-High Energy Cosmic Ray (UHECR) composi- 
tion at the highest energies, above E ~ 10 19 eV, is still 
a matter of debate. The measurements from the Pierre 
Auger Observatory (PAO) suggest a shift towards heav- 
ier nuclei [3]. The muon data from the Yakutsk EAS Ar- 
ray were also found to be consistent with the PAO mea- 
surements [4]. On the contrary, the measurements from 
HiRes [5] as well as the preliminary results of Telescope 
Array [6] are compatible with a proton composition. 

The cosmic ray spectrum shape [7, 8] is compatible with 
either a proton or an iron nuclei composition, for E > 
10 19 eV. Heavy nuclei may be accelerated to higher ener- 
gies than protons, which may lead to a transition towards 
heavy nuclei at the end of the spectrum. Additional in- 
formation on the primary composition may be drawn from 
the difference of propagation between protons and heavy 
nuclei in the Galactic Magnetic Field (GMF). Until now, 
most studies of UHECR propagation in the GMF dealt with 
proton or light nuclei primaries, with the exception of Ref- 
erences [9, 10, 11, 12, 13, 1,2]. 

In this work, we study the propagation of UHE iron nuclei 
with energies E > 60 EeV in the GMF. Both the impacts 
of its regular -large scale- and turbulent -small scale- com- 
ponents are discussed. Quantities in given directions of the 
sky, such as the deflection angles on the celestial sphere 



of UHE nuclei between their sources and their arrival di- 
rections at Earth, are strongly model dependent. However, 
quantities averaged over the whole sky, such as the mean 
deflection angles, are considerably less model dependent 
for all tested recent GMF models. In this sense, we are able 
to present generic conclusions. We take here the Prouza 
and Smida (PS) model [14, 15] for the regular GMF. For 
the turbulent component, we construct it as a superposition 
of a finite number of plane waves, following Refs. [16, 2]. 

Section 2 deals with the image formation of UHE nuclei 
sources in the GMF, and discusses the possibilities and 
challenges of "UHECR astronomy" in case a heavy compo- 
sition would be confirmed in the future. Section 3 presents 
a study of (de-)magnification effects of individual UHE 
iron nuclei source fluxes. Such effects are due to mag- 
netic lensing in the GMF, and can be considerable in a non- 
negligible fraction of the sky for heavy nuclei primaries. 

2 Images of UHE heavy nuclei sources 

The shape of images of UHE proton point sources on the 
celestial sphere has been discussed in detail in Refs. [17, 
18]. Except in directions to the Galactic plane and cen- 
ter, UHE protons are approximately deflected on the sky 
from their sources as 1/E. At the highest energies, outside 
the Galactic plane, there is one-to-one correspondence be- 
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Figure 1 : Image of the Centaurus galaxy cluster emitting 
UHE iron nuclei (from E — 60 to 140 EeV) deflected in the 
PS regular GMF model. Black disk for Centaurus position 
on the sky. Colours according to energies, see key. 

tween UHECR entering directions in the Galaxy and their 
directions on the sky as measured at Earth. 

For heavy nuclei primaries, there is no one-to-one corre- 
spondence in large regions of the sky, even at the highest 
energies. In such regions, a given extragalactic source can 
have several images on the celestial sphere. In practice, nu- 
clei from different images enter the Galaxy from the same 
direction, but at different positions in the physical space. 

Before entering our Galaxy, UHECR may be deflected sig- 
nificantly in extragalactic magnetic fields (EGMF) which 
are expected to correlate with the Large Scale Structure 
(LSS) of galaxies but are otherwise poorly known. If 
UHECR deflections in EGMF are as low as in the Dolag 
et al. scenario [19, 20], deflections of UHE iron nuclei 
would be negligible outside galaxy clusters. However, for 
a source located in a cluster, UHE nuclei would still be suf- 
ficiently deflected within the cluster to make it shine as a 
whole and make it appear as an extended source on the 
sky [21]. We compute in Fig. 1 the image of the Cen- 
taurus galaxy cluster in the PS regular GMF model, for 
energies E = 60 — 140 EeV, assuming the cluster con- 
tains one or several UHE iron nuclei source(s). The clus- 
ter is denoted by a black disk on the sky, and its cosmic 
ray arrival directions as measured at Earth are represented 
by small colour dots. If one also adds the contribution of 
the turbulent GMF, the UHECR arrival directions on the 
sky would be more spread. In some circumstances, its im- 
pact is less trivial than for UHE protons, especially at low 
Galactic latitudes [2]. Due to the poor knowledge of the 
GMF [22], Fig. 1 should not be considered as a predic- 
tion. The exact image of a given source is strongly de- 
pendent on the GMF model. Even the two most recent 
benchmark models proposed in Ref. [23] would result in 
different images. Nonetheless, the generic features visible 
in Fig. 1 are quite typical for UHE heavy nuclei sources, 
in all Galactic magnetic field models. In most regions of 
the sky, heavy nuclei sources have several images. Some 
of them may be strongly distorted and (dis-)appear above 
given energy thresholds, such as the image at high latitudes 
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Figure 2: Images of the supergalactic plane, regarded as a 
±10° band in supergalactic latitude &sg> emitting 60 EeV 
iron nuclei deflected in the GMF. PS model for the regu- 
lar GMF. Upper panel: No turbulent field added; Lower 
panel: Additional turbulent component with Bq = 4 /iG, 
L c = 50 pc and z$ = 3 kpc. The blue points denote the CR 
arrival directions. The grey lines stand for the supergalac- 
tic plane and the 6sg = ±10° circles. Map centered on the 
Galactic center, with increasing I from the left to the right. 

in Fig. 1. The energy ordering of events on the sky may 
be very far from the 1 jE ordering for protons, and in some 
cases, lower energy events may be closer to their sources 
than higher energy events. Images of more sources, and 
in other GMF models, are presented in Refs. [1, 24]. As 
shown in Refs. [24, 25], in some GMF models and for 
some sources, at least one of their images may still look 
like a proton image roughly enlarged by a factor Z, for nu- 
clei of charge Z. In such cases, one may still be able to 
detect them with the algorithm presented in Refs. [24, 25]. 
However, in most cases, a better knowledge of the GMF 
than currently available is needed so as to find efficient al- 
gorithms of source detection. 

If UHECR deflections in EGMF are as large as in the Sigl 
et al. scenario [26, 27], UHE heavy nuclei would still 
preferably be deflected within the LSS, with negligible de- 
flections in the voids of the LSS. Therefore, even in this 
scenario, the UHECR directions when entering our Galaxy 
should still correlate with the LSS, even if they would not 
point back towards their sources -or clusters. As a first ap- 
proximation, we study here the image of the supergalac- 
tic plane as a whole. Figure 2 (upper panel) presents the 
60 EeV iron nuclei image of the supergalactic plane, in the 
PS regular GMF only. Dark blue dots correspond to the 
arrival directions of cosmic rays on the sky. We consider 
the supergalactic plane as a ±10° strip in supergalactic lat- 
itude &sg- Grey lines stand for & SG = -10°, 0°, 10°. Fig- 
ure 2 (lower panel) presents the same image when a non- 
zero turbulent component is added to the PS regular GMF 
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model. We take here a turbulent field of correlation length 
L c = 50 pc. We assume a Kolmogorov spectrum with vari- 
ation scales between L m i n = 20 pc and L max = 200 pc, 
but for the chosen rigidities, our results do not depend on 
the magnetic field power spectrum. We take a spatial pro- 
file as in Ref. [18] with a rms strength at Earth B = 4 /iG 
and extension in the halo zq = 3 kpc. For such parame- 
ters, the UHE iron nuclei arrival directions are significantly 
more spread on the sky, making it more difficult to recog- 
nize the initial supergalactic plane image without turbulent 
GMF. Since the turbulent component is still poorly known, 
the dependence of results on its parameters have been in- 
vestigated in Ref. [2], 

3 Magnetic lensing in the GMF 

Magnetic lensing of UHECR in the GMF results in the (de) 
magnification of source fluxes as seen at Earth. The flux de- 
tected at Earth for a given source is changed by a factor A, 
compared to the flux one would have received if the GMF 
were set to zero. We call A the amplification factor. De- 
pending on the direction on the sky, it may be either larger 
(magnification) or smaller (demagnification) than 1. For 
UHE proton sources, this effect is negligible in all GMF 
models (A — 1), except for the unfavorable directions to- 
wards the Galactic plane and center. On the contrary, if 
UHECR are heavy nuclei, such effects would be substan- 
tial even at the highest energies E > 60 EeV. 

We compute A by backtracing a few times 10 5 iron nuclei 
from the Earth to outside the Galaxy. The densities of out- 
going nuclei are proportional to A -see References [1,2]. 

We present in Fig. 3 the logarithm of the amplification fac- 
tor, log 10 ("4)> for 60 EeV iron nuclei sources, depending 
on their positions on the sky. For the GMF, we take the PS 
model, with no turbulent component added (Bq = /iG). 
The black regions correspond to regions where log 10 (A) < 
—2. Let us call them "blind". In such directions of the sky, 
fluxes of 60 EeV iron nuclei sources would be demagnified 
more than 100 fold. Sources in blind regions cannot be de- 
tected neither by present nor near future experiments. In 
a few small regions (in white), the amplification factor is 
larger than 10. Only ~ 55% of the sky has flux modifica- 
tions by less than a factor ~ 3 (—0.5 < log 10 (.A) < 0.5). 

The exact amplification factor in a given direction of the 
sky is strongly model dependent. However, the distribution 
of regions with given amplification is less model depen- 
dent, for all recent GMF models. 

In Figure 4, we present the same plot as in Fig. 3, with 
a non-zero turbulent component added to the PS regular 
GMF. The turbulent field parameters that are considered in 
this figure are identical to those of Fig. 2 (lower panel). 
One can see that the regions of extreme magnification or 
demagnification (in bright and dark colours) tend to shrink. 

We compute the extension on the sky of regions with given 
amplification factors, following the method of Ref. [1]. 
Figure 5 presents the results in bins of width A log 10 (A) = 
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Figure 3: Logarithm of the amplification factor A of the 
flux of extragalactic 60 EeV iron nuclei sources, depending 
on their position on the sky. PS model for the GMF. Black 
regions correspond to log 10 (A) < —2. Plot in Galactic 
coordinates, with the Galactic anti-center in the center. 
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Figure 4: Same as in Fig. 3, with an additional turbulent 
component with Bq = 4 /iG, L c = 50 pc and zq = 3 kpc. 
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Figure 5: Fractions of the sky outside the Galaxy with am- 
plifications A. PS regular GMF model for all four his- 
tograms. Solid red line for no additional turbulent field. 
Green, blue and magenta lines for additional turbulent 
components with parameters respectively set to (Bq = 
1 /xG, zq = 3 kpc), (B — 4 /iG, z = 3 kpc) and (B = 
10 /iG, zq — 10 kpc). Correlation length set to L c = 50 pc. 
Bin of amplifications A below 10~ 2 for blind regions. 

0.5, both for the cases of Fig. 3 (no turbulent field - red solid 
line) and Fig. 4 (additional turbulent field with Bo = 4 /iG 
and zq = 3 kpc - blue thin dashed line). We also show 
the case of turbulent fields with (Bq = 1 /iG, zq = 3 kpc) 
and (Bq — 10 /iG, zq — 10 kpc), respectively with green 
thick dashed and magenta dotted lines. This plot confirms 
the above statement. For turbulent components that are 
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stronger and more extended in the halo, regions of strong 
(de-) magnification globally shrink. For instance, the size 
of blind regions is divided here by ~ 3 in the case of a 
turbulent field with B = 4/iG and z = 3kpc. A strong 
and extended turbulent component would however not nec- 
essarily be helpful for the search of UHE heavy nuclei 
sources due to the larger blurring of their images -see Sec- 
tion 2. For the extreme case (B — 10 fiG, z = lOkpc), 
the fraction of the sky satisfying —0.5 < log 10 (^4.) < 0.5 
reaches ~ 80%. The maximum and minimum amplifica- 
tion factors also respectively decrease and increase when 
the turbulent field strength and spatial extension are in- 
creased. 

4 Conclusions and perspectives 

We have investigated here UHE heavy nuclei propagation 
in models of the Galactic magnetic field, taking into ac- 
count both its regular and turbulent components. We have 
shown that several important effects appear in case of a 
heavy primary composition at the highest energies. 

We computed in Section 2 the images of an individual 
galaxy cluster and of the supergalactic plane, emitting UHE 
iron nuclei deflected in GMF models. We pointed out that 
for such a composition, non-trivial effects often appear. For 
instance, multiple and distorted images, or energy ordering 
of events far from the 1 /E behaviour expected for proton 
sources, can be seen for sources located in large fractions 
of the sky, in all recent GMF models. 

Except in a few favorable cases, a better knowledge of 
the GMF than currently available is still necessary in or- 
der to identify UHE heavy nuclei sources. Future radio 
telescopes, such as LOFAR and SKA, will enable one to 
improve the knowledge on the regular and turbulent GMF 
components, both in the Galactic disk and halo. 

In Section 3, we studied the (de-)magnification of UHE 
iron source fluxes depending on their positions on the sky. 
This effect is due to magnetic lensing of UHECR in the 
GMF. While it is negligible at the highest energies for pro- 
ton primaries, it becomes substantial for heavy nuclei. We 
showed that, at 60 EeV, sources located in up to about one 
fifth of the sky can have their fluxes demagnified by more 
than a factor 100 at Earth, making them undetectable by 
present and next generation UHECR experiments. This 
fraction is reduced for stronger turbulent components. For 
a turbulent field strength of 4 ^G at Earth and 3 kpc exten- 
sion in the halo, the fraction of invisible sky in the PS reg- 
ular GMF model is divided by more than three. This would 
however not facilitate source detection and reconstruction 
because of the larger blurring on the sky of source images 
due to the turbulent component. 
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